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We investigate the gravitational effect of a scalar field within scalar-tensor gravity as an alternative 
of the dark matter. Motivated by results of chameleon models, f(R) gravity and symmetron models, 
we study a phenomenological scenario where the scalar field has a mass and a coupling constant to 
the ordinary matter which scale with the local properties of the considered astrophysical system. We 
analyze the compatibility of this "alternative gravity" scenario at galaxy and galaxy cluster scales. 
Main results are: 1) the velocity dispersion of elliptical galaxies can be fit remarkably well by a scalar 
field, with model significance similar to the one obtained if a classical Navarro- Frenk- White dark halo 
profile is considered; 2) in particular, the analysis of the stellar dynamics and the gas equilibrium 
on elliptical galaxies has shown that the scalar field can couple with ordinary matter with different 
strength (different coupling constants) depending on the clustering state of matter components; 3) 
spiral galaxies and clusters of galaxies combined together show evident correlations among theory 
parameters (coupling constants and scalar field interaction length) which suggests the generality of 
the results at all scales and the way toward an unification of the theory for all gravitating systems; 
4) the gravitational effects of the scalar field and its viability as an alternative to the dark matter are 
confirmed by some preliminary test on strong lensing at galaxy cluster scales. 

Subject headings: dark matter - galaxies : elliptical and lenticular - galaxies : clusters : intracluster 
medium - galaxies : spiral - gravitation 



1. INTRODUCTION 

Dark matter and dark energy are nowadays widely ac- 
cepted as the main components of our Universe, although 
reliable clues about their origin, nature and properties 
are still missing. 

There is a long list of possible dark matter candidates 
ranging from standard to sterile neutrinos, from axions 
to super-symmetric c andidates, from light t o super- heavy 
scalar fields (see e.g. iBertone et al.l [20051 : iKraussI [20061 
and references therein) , while dark baryons seem to have 
a minor contribution. 

For what it concerns the nature of dark energy, the 
largest component in the mass-energy balance of the Uni- 
verse (« 74%), a coherent cosmological model explaining 
all open issues and the related observed phenomena has 
still to be found. 

These two "dark" ingredients are the pillars of 
the current cosmological "concordance" model, the 
ACDM. where A s t ands for the cosmological constant 
([Carroll et al.l 119921 : iSahni fc Starobinskj|2000T ) which is 
assumed to drive the Universe expansion, and CDM 



stands for Cold Dark Matter. This mo del provides a 
good fit to most of the cosmological data (Sanche z et al.l 
[2001 ISeliak et all [2001 iTegmark et all 120041 ) but it is 
also well known that it is affected by many serious 
theoretical problems that motivate the search of more 
general candidates for the acceleration-driver counter- 
part. Such models range from scalar fields rolling down 
self-interaction potentials to phantom fields, from phe- 
nomenological unified models of dark ener gy and dark 
matter to alternative theories of gravity (ICapo zzielk 
2002t iKoivisto fc Mota 1 120061: ILL Barrow fc Mota! 1 2007 
Capozziello fc Francaviglia 2008; Copcl and et alj |200C 
Manera fc Mota Il2006t i Mota Il2008t iPadmanabhanl 200c 



Peebles fc Rathrall2003 IKoivisto Mota fc PitrouH200S 
Capozziello fc De Laurentisll201lD. 

Among all these candidates, in IMota et al.l ((201 1L Pa- 
per I, hereafter) we have focused on a particular scenario, 
where a scalar field might be used to unify the cosmolog- 
ical scale acceleration of the Universe with the formation 
and dynamics of gravitational structures by mimicking 
dark matter on astrophysical scales. 

Scalar fields play an important role in connecting 
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cosm ology and particle physics (|Binetruvl 120061 : iLindd 
2008). In particular, we have s hown that theories 
like th e chameleon field th e ories (iKhourv fe Weltmanl 



[200l iMota fc Shawl 12006) IBrax et al.1 12004 
gravity (iHu fc Sawickl [2007 



f(R) 
models 



_ iStarobinsk 

iCapozziello et al.l 20031) an d symmetron 
(iHinterbichler fe Khomyf 120101: lOlive fc Pospelovl 120081: 
iDavis et al.l 120 111) , all share the fundamental ingredient 
of a scalar field which couples to matter via gravita- 
tional interaction, being short-ranged in highly dense 
regions, and long-ranged in low density regions. Such 
a scalar field would be relatively light on cosmological 
scales, thus describing cosmological evolution without 
differentiating much from the ACDM model, and, at the 
same time, it would also be able to satisfy local gravity 
constraints. Moreover, at different astrophysical scales, 
the effect of the scalar field would be suppressed or 
enhanced according to the local astrophysical density. 

Starting from these broad properties, we propose a new 
"paramctrization" for a massive scalar field theory where 
wc introduce a field mass (or equivalently an interaction 
length) and a coupling constant with ordinary matter 
which may change with scale. Our primary goal is to test 
whether this paramctrization, based on the chameleon 
or symmetron theories, can account for different obser- 
vational facts, regardless the underlying physical mech- 
anism producing such a scalar field. So me recent pa- 
pers (jOvaizu et al.l [20081 : IBrax et al.ll2012h analyze sym- 
metron and f(R) gravity theories (of which chameleon 
models are an extension) in the context of structure for- 
mation but they are considered as a background cosmol- 
ogy. Here we propose a scaling mechanism for which they 
can both explain dark energy (on cosmological scales) 
and dark matter (on astrophysical scales). 

In Paper I we studied the feasibility for the scale- 
dependent scalar field to work on different gravitational 
scales by using various cosmological indicators: type la 
supernovae and their Hubble diagram, low surface bright- 
ness spiral galaxies and their rotation curves, clusters of 
galaxies and their mass profiles. 

Here we extend our analysis to elliptical galaxies and 
model the velocity dispersion profiles. As we will show 
later, this will lead us to new results that have made 
necessary to re-analyze the gravitational objects studied 
in Paper I (spiral galaxies and clusters of galaxies) and 
finally unify all evidences in a more general scheme. 

The article is organized as follows: In Sj2]we give a brief 
but exhaustive summary of all the main properties of the 
scalar field theory and describe all the main hypothesis 
underlying our work. In Sj3] we accurately describe the 
used astrophysical data and the way we involved them in 
our analysis. In S|4] we show results concerning elliptical 
galaxies and in $5] we discussion the implications for a 
unified picture of all the gravitational structures we have 
considered. Conclusions are drawn in ST7] 



2. THE SCALAR-TENSOR FIELD THEORY 

The most general action g overning the dynamics of a 
scala r field can be written as (Esposit o-Farese fc PolarsO 
l200lh : 
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d^x^—g {F(cb)K - Z^){d$) 2 - V ((/>)} 



d 4 X C r 



(i) 



where g is the determinant of the metric g^ u , 1Z is the 

Ricci scalar, ipm are the various matter fields, C m is the 
Lagrangian density of ordinary matter, <p is the scalar 
field, F(tfi) and Z{(f>) are some functions of the field 
which regulate its dynamics, and V(4>) is the scalar field 
potential. Depending on the expression of the given 
functions F((f) and Z((f>) and of the potential V(<fi) 
one can recover "general" scalar field theories, like the 
chameleon or the symmetron mechanisms. The main 
consequence of this action is derivable when considering 
linear perturbations of matter and their related equa- 
tion: a wide class of theories leads to a perturbation 
equation like this (IDavis et al.l [20TTI: IBrax et al.l 120041 : 
iBertschinger fc Zukinll2008l ): 



G 



N 



(2) 



where of course the dots means time derivatives, k is the 
wavenumber length, a is the scale factor, H is the Hub- 
ble function, S m is the matter density contrast and Q m 
is the matter density parameter. The quantity G e // can 
be interpreted as an effective gravitational constant: in 
general relativity the gravitational coupling Gat is con- 
stant while in many alternative theories of gravity, the 
strength of gravity varies with time and also with place. 
A time-varying gravitational coupling is a well known 
properties of scalar-tensor theories and a generic fea- 
ture of all modified gravity theories where the Newtonian 
pote ntial and the spatial curvature potential are differ- 
ent (IBertschinger fc Zukinl 120081: lAcquaviva et al.1 120051: 
IClifton et al.1120051) . Within the context of massive scalar 
field models it has the general expression CGannoui i et al.l 
2009): 



G ef f(a;P,m;k) = G N [1 + 2/3 2 



1 



k 2 



(3) 



where /3 is the coupling constant of the scalar field with 
matter and, if the field is at the minimum, the scalar field 
mass is m 2 = V^. In particular the term proportional to 
to 2 results from the scalar field-mediated force, which is 
negligible if the physical length scale of the perturbation 
is much larger than the range of the scalar field-mediated 
force, namely, ifa/fc^m -1 . In this case the left hand 
side of Eq. ([2]) is well approximated by 3Q 2 n H 2 5 m /2 and 
the matter fluctuations grow as in general relativity. 

Taking the inverse Fourier transform of Eq. ([3]) it is 
straightforward to obtain the corresponding expression 
of the gravitational potential for a point mass distribu- 
tion, i/j(r). Remembering that a potential a i in real 

space yields a k~ 2 term in Fourier space, we can recog- 
nize in Eq. ([3]) the point-like gravitational potential per 
unit mass: 



^(r) = — (l + 2/3 2 e- mr ) = 

_G 
r 



(l + 2/3V 



r/L 



(4) 



where m is the mass of the scalar field, Lamm 
the interaction range of the modified gravitational po- 
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tential, i.e. the length where the scalar field is effective, 
and /3 still being the coupling constant between matter 
and the scalar field. The gravitational potential given in 
Eq. (j4|) has been derived for a point-like source, and it 
has to be generalize to extended systems in a numerical 
way. Depending on the particular gravitational system 
we consider, we will adopt different geometrical hypoth- 
esis: for cluster of galaxies and elliptical galaxies we will 
adopt the spherical symmetry, while spiral galaxies will 
be assumed as thin disks. 

The point-like potential can be split in two terms. The 
Newtonian component is: 



ipN(r) = - 



GM 



(5) 



and its extended integral is the well-known expression: 



*jv(r) 



GM(< r) 



(6) 



where M(< r) is the mass enclosed in a sphere with 
radius r. The correction term coming from the scalar 
field is: 

^(0 = ^(2^) , (7) 

so that to obtain the extended integral we have to calcu- 
late the following integral: 

*c(r) = / r' 2 dr' / sin(9W / du'il>c{r') , (8) 
Jo Jo Jo 

where the angular part is analytically derivable, while the 
radial integral has to be numerically estimated once the 
mass density is given. A fundamental difference between 
the corrected and the Newtonian terms is that in the lat- 
ter the matter outside the spherical shell of radius r does 
not contribute to the potential, while in the former the 
external mass distribution enters into the the potential 
integral, although the contribution is almost negligible in 
general. 

We also observe that a possible dependence of the cou- 
pling constant with scale, i.e. /3 — /3(r), must be consid- 
ered when evaluating the integral or the physical observ- 
able quantities that we will define in following sections. 
As we do not know what is the possible analytical behav- 
ior of /3(r), in our approach we will consider this quan- 
tity as a constant or weakly depending with the scale, 
i.e. d/3/dr 0, as we will verify a posteriori. 

2.1. Hypothesis 

In this section we want to describe in more details 
what are the main properties and requirements of our 
parametrization and the motivation of our approach and 
finally discuss the comparison with observations. 

In Paper I we performed our analysis by assuming 
one single scalar field working for any considered grav- 
itational scale. This scalar field was characterized by 
a mass (or interaction length L), which should be re- 
lated to the dimension of the gravitating structure under 
exam, and by a coupling constant /3 with the "kind of" 
matter which constitutes the object. We have started 
by the hypothesis that the matter was made only by 
the observed baryons (hot gas and galaxies in clusters of 
galaxies; gas and stars in spiral galaxies) with the scalar 



field generating a dynamical effect similar to the "classi- 
cal" dark matter. In practice, we replaced the eventual 
new and exotic dark matter component with an effective 
mass induced by modified gravitational interaction from 
the scalar field with ordinary matter. 

It is important to stress that the coupling constant /3 
was assumed to be unique and fixed for all the mass com- 
ponents. As we will discuss in jjU this turns out not be 
the case for elliptical galaxies, where we will consider the 
possibility that the scalar field coupling constant shall 
have different values depending on the different mass 
components of a galaxy. 

This eventuality might have two implications: 1) we 
have two separate values for the coupling constants each 
for any baryonic mass component (i.e. stars and gas), or 
2) we have only one coupling constant but its measure- 
ment might be affected by the matter state. 

It is possible to demonstrate that these two options are 
only apparently conflicting with each other. 

Indeed, in the first case the scalar field theory predicts 
that the field can couple i n different ways w ith different 
kinds of matter (see e.g., (jBrax et all 120041 ) ). In a clas- 
sical approach (i.e. general relativity plus dark matter) 
this would mean that the scalar field can couple differ- 
ently with ordinary matter (baryons, neutrinos, quarks, 
and so on) and with dark matter. 

In the second case, the different measurements depend- 
ing on the clustered states of the matter can be a conse- 
quence of a screening effect which can suppress the field 
effects and produce an apparent (measured) lower value 
of the coupling constant j3. Such a screening effect is 
called "thin-shell effect" in the chameleon theory and a 
similar effect is also present in the symmctron theory. 
It mainly affects gravitational systems where the inner 
value of the scalar field is different from the background 
one and also reflects in a compariso n between the in- 
ner and the external matter density (iBrax "etalll200l 
iCapozziello fc Tsuiikawa Il2008t iDavis et alj|20lij ). 

We suggest a possible mixed scenario to make co- 
existent the two options above. 

Let us take the case of a cluster of galaxies and the 
galaxies and the gas inside it (or star and gas in a galaxy 
on smaller scale): if the scalar field scales with the den- 
sity, we can argue that there are two scalar fields, one 
driving the formation and the dynamics of the cluster 
and another one driving the formation and the dynamics 
of the galaxies inside it. We can think that all the sub- 
structures inside the cluster experience the cluster-scale 
scalar field but, as long as the systems evolve, there is a 
point where the (over-) densities representing the galax- 
ies are as large as the one of the cluster itself and a 
screening effect arises and wash the cluster-scale scalar 
field effect out. This process would result in a suppres- 
sion of the coupling constant of the cluster-scale scalar 
field with matter in galaxies, such as these latter can be 
considered as clustered small structures in a larger one. 

In this picture, the more diffuse hot gas component 
mainly present at the cluster scale, can be considered as 
having only one typical scale, the one of the cluster. This 
assumption is valid provided that the hydrostatic equilib- 
rium is realized. In this case, the scalar field-gas coupling 
constant should contain information about the cluster- 
scale scalar field. The same assumption cannot be made 
for galaxy systems, as spiral and elliptical galaxies, where 
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the gas can be strongly disturbed by local phenomena 
(as stellar winds, supernovae, radio jets by active nuclei 
etc.). 

Within this framework, the total potential can be writ- 
ten in two equivalent ways. We can define: 



(9) 



if we want to separate the correction that the scalar field 
provide to the gravitational potential from the classical 
Newtonian term. Alternatively, we can write: 



*(r) = ^star(r) + * gQS (r) 



(10) 



where the suffix star refers to star component in galaxies, 
but it can be replaced by galaxy when writing the cluster 
potentials. Each of the two terms is finally given by: 



1W0) = ^N,star(r) + ^C,star(r; Pstar,L) , (f f ) 



and 



(12) 



where we point out that the same scalar field (with mass 
oc L^ 1 ) can interact differently with ordinary matter 
(two values for the coupling constant, p sta r and f3 gas ) 
depending on their clustering state. 

Finally, we have to remark here that our approach 
is implicitly making a static assumption for the gravi- 
tational structures we consider. This means that, for 
the moment, we are ignoring that: 1) the values of the 
scalar field parameters could be subjected to temporal 
evolution, so that an analysis of how perturbations and 
overdensity collapse work would be necessary; and 2) a 
dynamical analysis under the influence of a scalar field 
should be performed to verify the stability of such grav- 
itational systems. 

3. ANALYSIS PROCEDURES 

As pointed out in §1 , we want to test the scalar field 
hypothesis on galactic scale by starting with elliptical 
galaxies. 

Differently from spiral galaxies, whose easy inter- 
pretable flat rotation curves have been one of the first his- 
torical evidences for dark matter, elliptical galaxies are 
pressure supported systems dominated by hot random 
motions. The orbital distribution of stars is very diffi- 
cult to be modeled and consequently also the mass dis- 
tribution is highly uncertain, because of the well known 
mass-anisotropy degeneracy. 

One way to have some insight on the internal dynam- 
ics is to use the information stored in the line-of-sight 
velocity dispersion as a function of position by inverting 
the solving the Jeans equations. Under spherical symme- 
try, and assuming no-rotation, the only effective equation 
governing the galaxy equilibrium is the radial Jeans equa- 
tion: 



&{ta 2 r ) 
dr 



2^£a z r 
r 



d*(r) 



(13) 



where £(r) is the luminosity density of the galaxy, a r (r) is 
its radial velocity dispersion. The anisotropy parameter 
is defined as: 



Pa 



1 - 



where at is the one-dimensional tangential velocity dis- 
persion (defined as a combination of the two angular 
components of the velocity dispersion tensor a 2 = (a\ + 
a 2 ,)/2 and a r the radial component. When at = a r , the 
system is isotropic and /3 a = 0. Instead, when /3 a = I 
the system is fully radial anisotropic while for f3 a — > — oo 
it is fully tangential. 

In the Eq. (|I3j) . the unknown quantities are the 
anisotropy parameter and the mass which generates the 
potential, while £(r) is given by the tracer distribution. 
Thus, different combinations of orbital anisotropy and 
radial distribution of the mass can produce the same 
observed dispersion profile. This mass-anisotropy de- 
generacy can be solved by using independent measure- 
ment for the mass. One possibility for that is to use 
the information from X-ray emission from the hot gas 
(jMathews fe Brighenti Il2003l ). i.e. density and tempera- 
ture, to solve the hydrostatic balance within the galaxy 
potential, provided that the gas is at the hydrostatic 
equilibrium. This is not always true in elliptic al galaxies 
(jDiehl fe Statlerl [20071: iHumphrev et al.ll2006D . 

Despite all these modeling complications, Jeans 
analysis has been extensively used in elliptical 
taking great advantage of discrete kine- 
tracers probing the gravitational poten- 
t o many eff ec tive ra dii (R e ff). G l obular 
(iPuzia et all l200i iBergond et all 2006; 

~~ ti 12010 : 



galaxies, 
matical 
tial out 
clusters 



20091: IShen fe Gebhard 



Romanow skv et al _ _ _ .. 

Schuberthl l2010t iWoodlev et all 1201 01) or planetary 
nebulae (PNe , see e. g. INapolitano et all l200ll 120021: 



see e. g. UNap olitano et alJ \zwil \ 2UU2: 
[2001 iPeng et all 120041: iDougfas et all 



(14) 



iMendez et al.l . _ . _ _ 

20071: iCoccato et all I2009t iTeodorescu et all 1201 
Napolitano et al.l 120111 ) have made possible to extend 
kinematics up to 5-7 i? c ff. Furthermore, the analysis 
of satellite orbiting around th e galaxies could exten d 
up to 50-500 kpc (< 10 R eS ) (|Klvpin fe Pradal [20091) . 
However, while globular clusters have classically been 
used as mass tracers for bright galaxies, being their 
samples too small in ordinary elliptical galaxies, PNe 
have been systematically us ed to map the mass pro- 
file of ellipticals (see e .g. iRomanowskv et all 120031 : 
INapolitano et al.l 120091 INapolitano et al.l 120111 N+ll 
hereafter, and references there in). 

The selection of viable objects for our analysis has 
been made within the elliptical galaxy sa mple, observed 
with the Planetary Neb ula Spectrograph (iDouglas et all 
120021) and presented in (jCoccato et al.H2009t C+09 here- 
after) , where stellar kinematics of the central regions are 
combined with PN kinematics of the galaxy regions out- 
side the R c s- Adding PNe data has given strong hints 
about the mass profiles of elliptical galaxies, but has also 
opened much more questions and nowadays many alter- 
native sce narios are equally feasible . Tw o exemplary 
cases are (| Romanowskv et all [2003) and ( Dck el et all 
2005): in the former, using the PNe, the galaxies ve- 
locity dispersion profiles are found to decline with radius 
and dynamical modeling of the data indicates the pres- 
ence of little if any dark matter in these galaxies halos; in 
the latter, starting from disc-galaxy merger simulations, 
the lower than expected velocities are in fact compatible 
with galaxy formation in dark matter halos, thus depend- 
ing on inner dynamics (elongated orbits) or on projection 
effects. 
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For our analysis, we are interested on galaxies that 
have both extended (stell ar) kinematics and pu blished 
X-rays observations (e.g. iFukazawa et al.l [2006f ). This 
because we want to explore the dependence of the scalar 
field dynamics on the mass component of the gravita- 
tional systems under exam. 

From the 16 galaxies reported in C+09, the only galaxy 
which had a complete dataset for our purpose is NGC 
4374 since other galaxies for which both long slit and 
PN kinematics were available, like NGC 3377, NGC 3379 
and NGC 4494, did not have a reliable deprojected X-ray 
emitting gas density profil^B- All the needed photometric 
properties of NGC 4374 are reported in Table [1] 

3.1. Velocity dispersion analysis 
The general solution to the Jeans equation, Eq. (jT3j) , 



£(r)a 2 r (r) 



1 



f(r) J r 



' f(s)i(s)^ld S 



where the function / is the solution to 
din/ 



dlnr 



2fl,(r) 



(15) 



(16) 



By projecting the velocity ellipsoid along the line of sight 
one can derive the line-of-sight velocity dispersion, which 
is the kinematical quantity to compare to the observed 
equivalent ones reported in C+09: 



1 



I(R)af os (R) = 



R 2 



vr 



^=dr 

-R 2 



(17) 



ry/r 



R 2 



where R is the projected distance from the center of the 
galaxy and I(R) is the stellar surface brightness pro- 
file. To calculate the line-of-sight velocity dispersion one 
needs two ingredients: an analytical expression for the 
anisotropy function and the total gravitational potential. 
For what concerns the anisotropy function, the usual 
way of proceeding is to compare observations with pro- 
files derived from cosmological N-body simulations. One 
could use different forms: from the simplest hypothesis 
of isotropy (/3 a — 0), to a constant anisotropy profile; an- 
other commonly used model is the Osipkov-Merritt one, 
which however provides a poor fit to the simulations. 
We decided to w o rk with t he ani sotropy model given in 
IMamon fc Lokasl (12005A| |b1. 120061) : 



1 r 



2r 



(18) 



where r a is a typical anisotropy length, assumed to 
be r a ~ l4R e ff, this value thus providing a good fit 
to data from dissipation-less cosmological N-body sim- 
ulations (for a more exhaustive discussion about reli- 
able anisotropy models see Fig. 2 and Section 3.2 of 
(jMamon fc Lokasll2005Bl) . 

Then, for calculating the gravitational potential, we 
need to model the galaxy components, that in this case 
are: stars, hot gas and the central black hole. 

1 Fukazawa, private communication. 



The stellar luminosity density can be obtained by de- 
projecting the observed surface brightness profile; galax- 
ies in our sample are fitted with the well known and 
mostly used Sersic profile, i.e.: 



I(R) = Iq exp 



l/rri 



(19) 



where Iq is the central surface brightness (in units of L Q 
pc~ 2 ), a s the Sersic scale parameter and m the Sersic 
shape parameter. Starting from this relation, the lumi- 
nosity densit y can be obtained by the a pproximation first 
proposed in (|Prugniel fc Simie n f[l99l : 



£(r) = ht{r/a s ) , 



with: 



t{x) 



x v exp(— x 
Ltot 



l/m\ 



(20) 



(21) 



(22) 



47r m r[(3 — p)m]a i s ' 

where the function p i s defined in 
(jLima Neto. Gerbal fc Marqued [19991 ) as 

p~ 1.0- 0.6097/m+ 0.05463/m 2 , (23) 

and where the total galaxy luminosity (in solar units) in 
the V-band where observations for NGC4374 were per- 
formed is: 



4(B T -fi D -C B v- M B , q) 



(24) 



where Bt is the B-band apparent magnitude, Mq is 
the Sun absolute magnitude in the V-band, C'bv is the 
galaxy color, needed to convert all luminosity parameters 
from the band B to the band V (it is obtained from the 
extragalactic database HyperledtQ) and /irj is the modu- 
lus distance (given in Table [lj. We also need: 



Jo = 



2-7T m T\2m]a 2 



(25) 



which appears in Eq. (fT?]) through the surface brightness 
expression I(R). 

The central black hole is assumed to have a constant 
density inside its typical size, the Schwarzschild radius 
tbh i which is: 

2GnM B h ,„„n. 
Tbh = 5 i (26) 

where of course Gn is the Newton gravitational constant 
and c the light velocity. Then the black hole mass Mbh 
is obtained using the relation which exists between a su- 
per massive black ho le and the host galaxy luminosity 
(|Gultekin et "aT1[2009l) : 



M 



BH 



-^QS.gs+i.n-iogiicLtot/io 1 



(27) 



The density is assumed to be null outside tbh- 

The gas profil es are assumed to follow the trad itionally 
used /3-mode0 (Cavalie re fc Fusco-Femianolll978h : 



Pgas{r) 




3/V2 



(28) 



2 http://leda.univ-lyonl.fr/ 

3 The p g that appears here does not have any relation with the 



scalar field coupling constant. 
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TABLE 1 

Ellip tical galaxies. Column 1: Galaxy name. Co lumn 2: photometric b and. Column 3: modulus distance from ITonry et alI 
1200 111 and shifted by —0.16 mag as explained in ICoccato et alI 1120091) , Column 4: distance of galaxy derived from modulus 
distance. Columns 5: total B magnitude corrected for extinction and redshift. Columns 6: Sersic scale radius for a n = 4 

PROFILE (CLASSICAL De VAUCOULEURS PROFILE). COLUMN 7: SERSIC SCALE RADIUS DETERMINED FROM THE SERSIC FIT. COLUMN 8: 

Sersic shape parameter. Column 9: Stellar Surface brightness at a s . Column 10: maximum distance from the galaxy 

CENTER OF PNE DERIVED KINEMATICS (MAXIMUM DISTANCE FROM THE GALAXY CENTER OF PNE DETECTIONS). 



name band fjLjj D 


13 j~ 


4 

ai 


a 3 


m n s 


Rlast 


(Mpc) 




(kpc) 


(kpc) 


(mag arcsec - 2 ) 


(kpc) 


NGC4374 V 31.16 17.06 


10.01 


5.97 


9.34 


6.11 23.1 


25.23 (34.07) 



where the central gas density, p ga s,a, the gas core length, 
a g , and the const ant n are derived by fitting mass pro- 
files provided by (jFukazawa et al.ll2006t ). 

We remember that in the classical Newtonian approach 
the total dynamical mass is made of two counterparts, 
dark matter and baryons (stars, gas and black hole), so 
that: 



M tot (r): 



d* 



G dr 

rf_ ( d^ Nybar 
G { 



d* 



N.DM 



dr 



dr 



= M bar (r) + M DM (r) , 
where of course 'J at is the Newtonian potential and: 



(29) 



^N,DM PDM 
^N,bar Pbar ~ Pstar 



(30) 



Pgas + PBH 



As we pointed out in the § d2J), in our approach we have 
that the total gravitational potential is made of a New- 
tonian term and a corrective one: 



so that the dynamical mass is: 

Mtot{r) =G— = 

_ r 2 /d*Ar : d*c' 
~~ ~G (~dT 



(31) 



G V dr dr J 

= M bar (r)+M eff (r) 



(32) 



where the effective mass is due to the modification of 
gravity produced by the scalar field, instead of requiring 
a new kind of matter as the dark one. Of course, the 
term it c has contributions only from visible baryonic 
mass, i.e. stars, gas and black hole. 
By comparing Eq. ([29| and Eq. (|32|) it is straightforward 
to observe that if we want that the scalar field fits data 
as well as dark matter we need: 



M eff (r) ~ M DM (r) 



(33) 



For this reason we also realize a fit of our data in 
the classical context of dark matter in order to perform 
a comparison between the two approaches. We used 
the classi cal Navarro-Frenk- White (NFW ) model density 
given in ([Navarro. Frenk fc Whitdll996l ) by the relation: 



Pdm (r) = p^ 1 - 
\a>d 



r 

1 + — 

ad 



(34) 



In |M anion & Lokas 2005B) more dark matter models 
are given, one being the gen eralized NFW model intro- 
duced by (Jin g fc Sutoll2000l ). with inner slope —3/2 in- 
stead of — 1 as it is in the classical NFW profile; and the 
other one being the convergent model of (jNavarro et al.l 
I2004D , with an inner slope that is a power- law function of 
radius. We have verified that using these two models does 
not give any substantial change in the general mass pro- 
files and in the fitting of velocity dispersion curves, being 
the classical NFW model sufficient for our requirements. 
Generally, one is used to convert the quantities appear- 
ing in Eq. (jMf. i.e. (Pq M , ad) in more useful quanti- 
ties as the virial radius, r v , namely the radius enclosing 
a mass whose mean density is w 100 times the critical 
density of the Universe, and the concentration parame- 
ter, C v i r — — . 

Finally, with all these ingredients, we are able to de- 
rive the observed quantity, i.e. the line-of-sight velocity 
dispersion, ai os . We underline here that the data we 
consider are derived from the sum in quadrature of two 
terms: the line-of-sight velocity dispersion and of the 
rotation velocity. The final quantity, y af os + v 2 (see 
N+ll), is a more efficient indicator of the total kinetic 
energy and is also « oi os because of the rotation velocity 
in NGC4374 is not dynamically significant compared to 
random motion, being ~ 50 km/s -1 against a velocity 
dispersion of ~ 200 — 250 km/s -1 . 

Depending on the approach we consider, o\_ os will be 
a function of different parameters: in the classical ap- 
proach with a NFW density profile for the dark matter 
component it will be cr/ os = ai os (R; p® M ,0,4, Y*), while 
in the modified gravity approach by using a scalar field 
it will be oios — <Ji os {R;j3,L,Y*). The quantity Y* is 
the stellar mass-to-light ratio we need to convert stellar 
luminosity density to mass density. For the statistical 
analysis we will use the chi-square function defined as: 



Xell 



AT 

E 

.7=1 



{vios,th{Ri] {9i}) @los.obs 



(35) 



where N is the number of data, a 2 the measurement 
variances and {9f\ is the parameters theory vector which 
is, respectively, = {Pq AI , a<z, Y*} in the dark mat- 

ter approach, and {6i} — {(3, L, F*} in the scalar field 
approach. For minimizing the \ 2 we use the Markov 
Chains Monte Carlo Method (MCMC) an d test their con- 
verge nce with the method described by (jDunklev et al.l 
120051) . 

4. ELLIPTICAL GALAXIES: RESULTS 
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In this section we apply the scalar field approach to 
the elliptical galaxies. 

As discussed in Sj3j the equations of our dynamical ap- 
proach in the alternative scalar field scenario depend on 
the parameters that characterize the scalar field proper- 
ties, i.e. its interaction length L and the coupling con- 
stant /3 (or, as we will discuss, various values for it de- 
pending on the number of mass components of each ob- 
ject) and the stellar mass-to-light ratio, Y*, which we 
need for the conversion of the stellar surface photometry 
to the stellar mass density. 

Even if the used MCMC method makes it possible to 
establish some priors on the fitting parameters, we de- 
cided to leave them as free as possible; for this reason we 
have f3 > (as in all the expression it enter as (3 2 , we 
do not have the possibility to really distinguish between 
a positive or negative values; moreover, the scalar field 
theory predicts it to be positive), and L > and Y* > 
as their are positive definite quantities. 

To check the validity of our analysis and the ability of 
deriving information about the scalar field from the ob- 
served velocity dispersion as the only kinematical quan- 
tity to model, we start by performing an ideal experiment 
using a mock galaxy with all intrinsi c quantities fixed 
followi ng the same prescriptions of ()Mamon fc Lokasl 
I2005BD . 

According to their prescriptions we have considered an 
elliptical galaxy with a baryon and a dark matter com- 
ponent characterized by the following parameters: 

• B-band luminosity L*^ = 1.88 x 10 10 h^ Q 2 L Q 
(where h 70 = 1 if H ' = 100 • h = 70 km s" 1 
Mpc -1 ), from w hich, using all the relatio ns given 
in Section 2 of ([Mamon fc Lokasl [2005Bh . we ob- 
tained a Sersic shape parameter m = 3.12 and a 
Sersic length a s = 11.6 h~~ 2 pc; 

• a typical stellar mass-to-light ratio Y* = 6.5; 

• a total mass-to-light ratio Y = 100 corresponding 
to a virial radius, r V i r = 79 R e ff and to a concen- 
tration parameter c = 9.70; 

• a black hole to stellar mass ratio Mbh/M* = 
0.0015; 

• a /3-model for the gas component with index f3 g = 
— 1.5 and core radius r c — R e ff/10. 

The velocity dispersion of such mock galaxy is shown 
in the left panel of Fig. Q] as a dark dashed line. This has 
a slightly decreasing trend with radius which eventually 
change its slope outside 50-60 kpc, i.e. well beyond typ- 
ical radial coverage by PN kinematics. 

In fact, we tried to recover (only in a qualitatively 
way) the total velocity dispersion profile with a scalar 
field, using Eq. (fTT|) as shown in the left panel of Fig. [TJ 
The scalar field prediction with a single coupling con- 
stant ((3 = 0.05 and L « 1000 kpc) (dark dot-dashed 
line) for all mass components does not match the velocity 
dispersion profile of the dark matter case (dark dashed 
line). We can only use such values because the dispersion 
velocity in the inner region is completely dominated by 
the stellar component; it works like a sort of "normal- 
ization" factor and strongly constrains the value that ft 



can have. A value too much different and higher than 
this one, would give a completely wrong velocity disper- 
sion reconstruction. In particular, the scalar field total 
velocity dispersion is almost equivalent to the stellar ve- 
locity dispersion (in a classical dark matter scenario) well 
beyond the radial extent of PNe measurements (vertical 
dashed line), eventually rising only after this limit. 

On the other side, if one assumes a different coupling 
constant for the two mass components (stars and gas; 
flstar = 0.05, figas = 5.6 and L ss 90) kpc) it is pos- 
sible to decouple the effects from any of them; in fact, 
the model which includes two coupling constants (dark 
solid line) nicely reproduces the velocity dispersion pro- 
files of the classical case, i.e. the scalar field can mimic 
the profile of a NFW dark halo. We also specify in this 
context that changing the length parameter r a which ap- 
pears in the anisotropy function f3 a (r) from the chosen 
be st value, r a ~ 14R e ff, to r a ~ lAR e ff (as discussed 
m (M amon fc Lo kas 2005B|), beyond of spanning a wide 
range of values, does not change our results. 

The circular velocity curve, calculated from the rela- 
tion v 2 = r d^/dr, is instead rather flat out to 100 
kpc, clearly showing how dominant is the dark matter 
in the total galaxy potential. This latter property is 
what we expect to reproduce with the scalar field. We 
finally stress that if we had adopted the circular veloc- 
ity as an observational quantity (as it is done in spiral 
galaxies), we would not have any chance of discriminat- 
ing between the two approaches, i.e., one or two coupling 
constants for the scalar field and the baryonic matter. 
This is shown in the right panel of Fig. [T] where we can 
see how the two different cases for the scalar field give two 
equivalent reproductions of the NFW profile at least in 
the depicted distance range (the two representing lines, 
the solid and the dot-dashed one, are indistinguishable 
as perfectly overlapping) ; this is an important issue to be 
considered for when evaluating results for spiral galaxies 
in next sections. 

Starting from this preliminary considerations, we can 
move to the analysis of a real system: NGC 4374 (see 
Table [J). 

4.1. NGC4374 

We will make use of the PN sample discussed in N+ll 
which we refer the reader to for the details of the PN 
sample properties and the derivation of the kinemati- 
cal profiles we will use in our analysis; the PN disper- 
sion profile extends out to ~ 5Pe// which is a distance 
large enough to explore any deviation of the galaxy dy- 
namics from a pure Newtonian no dark matter behavior. 
N+ll has provided a full Jeans analysis of the system 
in the ACDM framework and showed that the dark mat- 
ter distribution in this galaxy is nicely consistent with 
a "standard" NFW halo. In this work, the authors use 
dispersion-kurtosis constraints in order to alleviate the 
mass-anisotropy degeneracy. 

This galaxy has been also analyzed in the context of 
f(R) theories (iNapolitano et al1l2012[ ). where a Yukawa- 
like modification of the gravitational potential is adopted 
as alternative to dark matter. This approach is different 
from the one we are adopting in the context of the scalar 
field where the potential allows us to break the contri- 
bution of all mass components (stars and gas) which is 
a peculiar feature of our theoretical scenario that cannot 
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Fig. 1. — Mock elliptical galaxy analysis. Left panel: velocity dispersion profile; Right panel: circular velocity profile. Dark lines: solid - 
scalar field with two coupling constants; dot-dashed - scalar field with one coupling constant; dashed - classical NF W dark matter profile. 
Light lines: dashed - stellar velocity dispersion (circular velocity) in the classical dark matter approach; dotted - gas velocity dispersion 
(circular velocity) in the classical dark matter approach; dot-dashed - dark matter velocity dispersion (circular velocity). Vertical lines: 
dotted - effective radius R e //, limit achievable with only stellar photometry; dashed - 5 — 6 R e //> limit achievable with PNe reconstructed 
kinematics; solid: virial radius r v i T . 



be included in their physical model. 

Before we proceed with the scalar field models as in 
$3j we have to underline a preliminary but important 
dynamical choice: the stellar component is the one that 
dominates the velocity dispersion profile, because it dom- 
inates the profile in the inner region and this implies a 
stronger weight on the global fit to the velocity disper- 
sion. For NGC 4374 there are different literature mod- 
els of its stellar photometry (mainly depending on the 
extension of the adopted datasets): a typical De Vau- 
couleurs profile ( a Sersic profile with ind ex m = 4), with 
R e ff — 5.97 kpc (jCappellari et al.l l 2006fk a Sersic profile 
with Re ff = 11.69 kpc and m = 7.98 (jKormendv et al.1 
2009); and a Sersic profile with R e ff = 9.34 kpc and 
m = 6.11 (N+ll). The first model has some problem 
in fitting data at low distances from the center and can 
also affect large radii data, as the De Vaucouleurs profile 
is not always able to describe all the intrinsic features of 
an elliptical galaxy. The second model allows a better fit 
to the stellar profile but fails to recover the behavior of 
surface brightness at small radii (R ;$ 0.5 kpc). Finally, 
the third model gives a very good fit of the stellar pro- 
file and in a wider range than the previous one, namely, 
0.09 < R < 38 kpc, out to the distances covered by 
PNe observations. For this reason we decided to adopt 
the model profile from N+ll in the following dynamical 
analysis. Here we adopted two different approaches in 
order to optimize at the best the goodness of the surface 
brightness reconstruction: a first one considering all the 
data, and a second one considering only data points with 
R > 0.09 kpc. 

We have started by using our model machinery assum- 
ing the standard Newton dynamics and a NFW dark 
halo and cross-check our results with the ones presented 
in N+ll. The results are in Table [2] which somehow 
differ with the results discussed in their § 3.5.1. Before 
we go into the details of the discrepancy, we need to re- 
mark that there are some critical differences between the 
two approaches. Here we use some parametrized stel- 
lar surface density, while in N+ll they use an interpo- 
lated function. Furthermore, we are assuming here an 



anisotropy profile which is rather different from the one 
constrained by N+ll where they have used the kurto- 
sis information. These two main differences can produce 
some substantial divergences in the modeling of the very 
central data points thus affecting the parameters which 
are more sensitive to the small radii fit. 

For this reason we have repeated the models including 
all the kinematics data point, and excluding the data 
at R > 0.09 kpc and R > 1 kpc and find substantial 
differences as shown in Table [2]) . 

In particular, including the whole sample we have ob- 
tained a lower value of the central density p® M and 
an higher value for the NFW radius with no statis- 
tically significant difference between the total and the 
cut sample, but if we consider the expected relation be- 
tween these two parameters derivable from a collisionless 
ACDM universe with WMAP5 parameters (Eq. (13) and 
blue contours of Fig. 7 in N+ll), we have correspondence 
at 1(7 level. The concentration parameter, c 4, matches 
with the lower limit of most of the cases shown in the Ta- 
ble 2 of N+ll, and in particular with their assumed best 
reference model (i.e., an adiabatically contracted NFW 
profile with an anisotropy distribution different from the 
one adopted here, which should result in a more con- 
centrated dark matter density profile). The virial radius 
is notably larger; while the virial mass is perfectly con- 
sistent with their results even if showing a wider confi- 
dence level extending preferentially to higher values. If 
we compare these results with the Fig. 11 in N+ll, we see 
that our results for the parameters (c^,-, log M„i r ) fall in 
the region limited by the results inferred from late-type 
galaxies dynamics a nd from weak lensing o f all type of 
galaxies and groups (jNapolitano et al.ll2009f ) ; in particu- 
lar the values coming from the r > 0.09 kpc sub-sample 
perfectly match with this last curve. For what it concerns 
the luminous matter parameter, we found values slightly 
higher for the stellar mass-to-light ratio, F* ~ 6.6, that 
is however perfectly compatible with a Salpeter Initial 
mass function (Fig. 5 in N+ll). 

Looking to the left panel of Fig. [2]we see how the stellar 
component dominates the velocity dispersion profile only 
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TABLE 2 

Elliptical galaxies: NFW dark matter. Column 1: x 2 type. Column 2: central NFW density. Column 3: NFW radius. Columns 
4: stellar mass-to-light ratio in the observation photometric band. Column 5: NFW concentration parameter. Columns 6: 
virial radius. Column 7: virial mass. Column 8: stellar mass-to-light ratio in the B band. 
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TABLE 3 

Elliptical galaxies: Scalar field. Column 1: \ 2 type. Column 2: Coupling constant of scalar, field and star component. Column 
3: Coupling constant of scalar field and gas component. Column 4: Coupling constant of scalar field and black hole 
component. Columns 5: Scalar field interaction length. Column 6: stellar mass-to-light ratio in the related band. 
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Fig. 2. — Left panel. Velocity dispersion curve of NGC4374. Dark lines: solid line is the total velocity dispersion in the scalar field 
approach, i.e. scalar field + baryons; dot-dashed line is the total velocity dispersion in the classical approach, i.e. dark matter + baryons. 
Light lines: dashed line is the stellar velocity dispersion in the NFW approach; dotted line is the dark matter velocity dispersion. Right 
panel. The same as before but in logarithmic scale. 



in the very central regions (r < 1.0 kpc), with the black 
hole contribution being important only for very small 
scales, while the NFW dark matter becomes dominant 
in the line-of-sight velocity dispersion profile at w 10 kpc 
which corresponds to w 1 - 2 R e / / ■ 

Based on this, when using only data with r > 1.0 kpc 
we have been able to recover results more similar to N+ll 
for what concerns the NFW profile: an higher central 
density parameter for the dark matter profile (even if it 
is again lower than the one measured in N+ll); a smaller 
value for the NFW length (but still higher than N+ll); 
the couple (c V i r , log M V i r ) is now compatible both with 
weak lensing inferred trend and with the relation derived 
from WMAP5-based simulations; finally, the virial radius 
is now only a 30% higher than their value. 

If we now move to the scalar field alternative approach, 
it is clear (right panel of Fig. [5]) that this is as success- 
ful as the classical dark matter approach in modeling 
the dispersion profile. The largest differences with re- 



spect the NFW profile are found both at the very small 
radii, where the scalar field model shows a steeper slope 
while the NFW one seems to reach a plateau, and at 
the very outer region, where the scalar field model stays 
flatter than the NFW profile, although more extended 
data would allow to adjust the two model better at the 
largest distance and possibly recover a better agreement 
also at shorter scales. However, as seen by the \ 2 re- 
sults in Table [31 the best fit is also in this case very good 
and while the NFW and the scalar field approaches are 
almost equivalent when using all the data points, if we 
considering the value of the reduced x 2 (with the NFW 
X 2 being slightly smaller than the scalar field one) for 
the best-fit to the data points with R > 0.09, the scalar 
field turns out to provide a far better significance of the 
fit with respect the NFW model. 

In this case, the stellar mass-to-light ratio is smaller 
than the NFW-based one and it is now more compati- 
ble with the Kroupa IMF (|Kroupall2001[ ) values found in 
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N+ll. Even more importantly, the scalar field parame- 
ters turned out to be consistent with what we argued for 
the double coupling constant hypothesis. 

First, the stellar mass component shows the lowest 
coupling constant among all the mass components, even 
lower than the black hole one. As anticipated, this can be 
the consequence of the "average" effect of the screening 
action made by the scalar field on the stellar component. 
The same does not happen to the black hole, since this 
represents a singularity and it is difficult to detail the 
change in the field from inside to outside and the com- 
parison with its classical Newtonian force. 

Second, the coupling constant which refers to the gas 
seems to be mostly correlated with the galaxy gravita- 
tional potential and the scalar field mass. 

To conclude, we stress that both the classical dark mat- 
ter approach and the scalar field seem to be unable to de- 
scribe the small "shoulder" in the dispersion profile that 
is present at R rj lOkpc, precisely where the PNe data 
overlap the only stellar kinematics. One may think that 
this can be a consequence of some unaccounted orbital 
anisotropy. In fact, this mainly depends on the choice we 
made to use a parametric form of the star density £(r) 
instead of a smoothed light profile like in N+ll where it 
was possible to recover all the details of the dispersion 
profile (with a small degree of radial anisotropy). This 
means that the observed kinematics is somehow strongly 
sensitive to the tracer space density distribution. 

5. UNIFIED PICTURE 

The result obtained with the elliptical galaxy NGC 
4374 is not meant to have the sufficient generality to 
drive any conclusion on the new proposed multi-coupling 
constant scenario. However, it is not surprising that the 
proper handling of physically more complex, multicom- 
ponent systems like ellipticals has opened new perspec- 
tives also in the scalar field approach. If our scenario 
is correct, this shall be in fact made more evident on 
gravitational systems where the evolutionary status is 
different (it is the case of spiral galaxies, as we will show 
in next sections) or where matter shows very different 
phases. But a more detailed analysis in this sense is 
out of the purpose of our work and will be developed in 
forthcoming papers. 

Since we do not expect to extend the elliptical sample 
any soon, due to the difficulty in finding galaxies having 
both extended stellar (PNe) kinematics and not faint gas 
profiles, in this section we try to re-interpret the results 
of Paper I in terms of the new evidences found for NGC 
4374, i.e. the possibility of disentangling the various mass 
components with respect of the coupling constant with 
the scalar field. In this respect we aims at drawing a 
common scenario for all the gravitating systems (galax- 
ies: spirals and ellipticals, and clusters) in a consistent 
way. 

5.1. Clusters of Galaxies 

By comparing Figs. [3] - [4] with the corresponding fig- 
ures in Paper I (Figs. 3 - 4 in that paper), we see that the 
the new approach of using the two mass component has 
allowed us to improve the model fit in 9 cases, especially 
at small radii. In one case (RXJ1159) the new approach 
has produced a remarkable fit where the model in Paper 



I failed completely. Only for 3 clusters the mass profile 
reconstruction has not got any benefit by the two com- 
ponent approach, although for all sample the deviation 
of the best fit from the data in never worse than at la 
confidence level. 

Looking at Table @] we can quantify better the differ- 
ences with Paper I: the interaction length slightly dimin- 
ishes in the most of the objects, still staying compatible 
with typical cluster scales (ra 1 Mpc; see also Eq. (36)), 
while the coupling constant with gas is generally larger 
than previous estimates. Two objects in particular have 
very large coupling constant values, i.e. clusters A262 
and MKW4. 

On the other side, the coupling constant related to 
galaxies is always much lower than the one of the gas 
one (as it is expected from the discussion in § 12.11 except 
for two cases, i.e. clusters A133 and RXJ1159. It is in- 
teresting to note that among all these exceptions, two of 
them (MKW4 and RXJ1159) are considered more simi- 
lar to group than to real cl usters of galaxies or e ven as 
extended elliptical galaxies (Vikhlini n et al.l 12005). thus 
for all these systems the argument adopted for normal 
galaxies might not apply straightforwardly. 

For the other two objects no particular features where 
found in literature to explain the anomalies. We can only 
remark that A262 seems to be as small as MKW4 (i.e. 
they share the same data extension) and differently from 
this latter is believed to be a normal cluster of galax- 
ies in literature. A possible relation between the clus- 
ter scale and the gas coupling constant is possible: if 
we think about this last one as a concentration param- 
eter, we can expect that smaller structures exhibit an 
higher value for it. For A133 nothing peculiar was found 
in literature so to justify the high value for the galaxy 
coupling constant we found; we can only verify that the 
fit with data is really good. As we did in Paper I we 
want to find some relations among the scalar field pa- 
rameters and the physical properties of the considered 
gravitational systems to establish if this alternative sce- 
nario can be a valid alternative to general relativity and 
dark matter. First of all, in the right panel of Fig. [5] 
we have a relation between the interaction length, L, 
and the radius rsoo, the distance from the center corre- 
sponding to an overdensity ps 500 times relative to the 
critical density at the cluster redshift. We prefer this 
quantity with respect of the virial radius r V i r used in Pa- 
per I because while the latter is deriv ed using a relation 
coming from cosmqlogica l simulations (| Bryan fc Nor maul 
Il998t lEvrard et al.lll996T) and thus depending on the gas- 
weighted aver age temperature of th e cluster, the former 
is derived in ([Vikhlinin et al.l [20051) directly from obser- 
vational data by using the hydrostatic equilibrium equa- 
tion. We have also verified that the two distances are 
proportional, having rsoo cx r^° 39 . 

The relation between the scalar field length L and the 
radius r5oo is derived from an error weighted fit exclud- 
ing the previously described four peculiar clusters (right 
panel of Fig. [5]): 

logL= (-1.59 ±0.82) + (1.49 ±0.27)- log r 500 ■ (36) 

We have also found a relation between L and the average 
gas- weighted temperature < T >, left panel of Fig. [5] 

logL = (2.45 + 0.11) + (0.68 + 0.17) -log < T > . (37) 
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Fig. 3. — Dark matter profile vs radii for clusters of galaxies. Dashed line is the observationally derived estimation of dark matter; solid 
line is the theoretical estimation for the effective dark matter component; dot-dashed lines are the 1-cr confidence levels given by errors on 
fitting parameters plus statistical errors on mass profiles. 



As we pointed out in Paper I, there exists a relation 
between the cluster mass and average temperature, i.e.: 

M A /T 3 / 2 oc H /H(z) , (38) 

where A is the overdensity level relative to the critical 
density at the cluster redshift, so that Migo = M V i r . 
With the previous phenomenological expressions holding, 
this mass-temperature relation can be easily converted 
in: 

L - 27 <xH{z)/H . (39) 

In Fig. |6] we can verify that assuming a fiducial WMAP 
quintessence mode0, with fl m = 0.259 and w = —1.12, 
and assuming H = 72 A km s _1 Mpc -1 , the previous 



relation is able to match the chosen fiducial model fol- 
lowing the relation H(z)/H = (0.169 ± 0.006)L 27 . 

In Figs. [7] we also represent the scaled temperature 
profiles versus the distance from the center of any clus- 
ter scaled with respect to the scalar field length obtained 
by fit, showing the same good reproduction of the self- 
similarity that characterizes the classical dark matter ap- 
proach and, as we found in Paper I, the absence of the 
subgroups which clusters ar e divided in depending on the 
mean temperature values (|Vikhlinin et al.l [2005) . The 
only exception with respect of Paper I is A133, which 
now appears to be to more stretched down, even if no 
peculiar elements have been found in literature and no 
problems in performing its analysis have been encoun- 
tered. 



4 http:/ /lambda. gsfc.nasa.gov/product/map/current/ 
parameters.cfm 



5.2. LSB spiral galaxies 
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Fig. 6. — Comparison between the Hubble function H (z) calculated from the fiducial model described in the text and the Hubble function 
obtained by the empirical relation H(z) <x L 0,27 . Error bars are calculated from the errors on the interaction length L. 




r/L 



1.4 



1.2 



1.0 



0.8 



0.6 



0.4 

























1 A \ 



0.01 0.02 



0.05 0. 10 0.20 
r/L 



0.50 1.00 2.00 



Fig. 7. — Temperature profiles for all clusters plotted as a function of distance from the center and in units of the scalar field length 
L. The temperatures arc scaled to the mean (gas-density weighted) cluster temperature. Dashed line is for high temperature clusters 
(< T > > 5 keV); solid line is for intermediate temperature clusters (2.5 < < T > < 5 kcV); dotted line is for low temperature clusters 
(< T >< 2.5 keV). 
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We now consider the case of spiral galaxies. The mod- 
eling of these systems involves the same number of pa- 
rameters of elliptical galaxies, i.e. the parameters associ- 
ated to the scalar field and the stellar mass-to- light ratio, 
Y*. As discussed for the elliptical case, the latter param- 
eter cannot vary arbitrarily but has to be consistent with 
stellar population models. As for spiral systems one can 
expect Y* to vary in the range « 0.5 and « 2 for the 
galaxy sample that we have adopted in this work, we 
have decided to be conservative and leave the model pa- 
rameters to vary in the interval [0; 5] (for any further and 
detailed description see Paper I). 

Before we proceed to illustrate the results, we need to 
remark that the two component approach, for the spiral 
galaxies, can suffer of some additional noise source with 
respect the analysis performed on paper I and can some- 
how affect the results. As we have exhaustively discussed 
in Paper I, the gas data are very noisy and in some cases 
they can also show particular features which are unavoid- 
able due to intrinsic dynamical properties (e.g. negative 
velocities due to counter rotating discs, large scatter from 
non-circular motion, etc.). While these effects have been 
mitigated in Paper I, where using stars and gas together 
allowed the more circular velocity of stars to dilute the 
noisy features of the gas, here the gas component to the 
scalar field might result highly uncertain and cause an 
higher value for the gas coupling constant than the one 
of stars. In order to improve the accuracy of the fit for 
the gas component, we have tried in this work to obtain 
a more accurate fit of the gas density: the result of this 
new analysis is evident in the smoother trend of most of 
the galaxies. 

The new models of the spiral rotation curves seems to 
give results that are less clear than the one obtained for 
clusters in the previous section (see Figs. [8] - [9j): by de- 
coupling the gas and the star components, we do not gain 
much in the quality of the fits (except a slightly better 
accordance to data of NGC4173 and NGC9211, where 
the models are closer to the ones obtained in Paper I). 
This is the consequence of the evidence showed in the 
previous section when testing the method on elliptical 
galaxies. We have seen that the effect of the decoupling 
of the mass components does not show up in the cir- 
cular velocity but only the velocity dispersion. In the 
LSB analysis, as we are modeling the circular velocity 
we should expect a minimal improvement by the addi- 
tion of a further coupling constant which implies that we 
are limited to reliably decouple the different mass com- 
ponents in spiral galaxies because of the unsuitability of 
the dynamical probe (the rotation curve). 

Of course, we remind that this inducted idea has to be 
considered with caution, since elliptical and spiral galax- 
ies are completely different gravitational structures from 
a dynamical point of view: elliptical galaxies are domi- 
nated by random motions while spiral galaxies are domi- 
nated by ordered (rotational) motions. Moreover we can 
add another element which make our global results more 
homogeneous and consistent each other: the decoupling 
(with an effective improvement in fits) is given for el- 
liptical galaxies and clusters of galaxies, which are both 
hot systems (i.e. dominated by random motions) and are 
supposed to be virialized, while this does not happen for 
spiral galaxies. Thus it is possible that the possibility 
of decoupling the different mass components is related 



to the reached equilibrium in those two kinds of gravi- 
tational structures, while in spiral galaxies stars and gas 
are again more strictly correlated. In that case, we would 
have also a correlation of our scalar field analysis with 
the evolutionary state of the gravitational objects; but 
to verify this is out of the purpose of this paper. 

For all these reasons we have proceeded with the star 
and gas component decoupling approach also for LSB 
systems and we have found we have a clear distinct be- 
havior for the coupling constants. First of all, unsurpris- 
ingly the gas coupling constant seems to be very similar 
to the only one coupling constant of Paper I, and it is 
very well constrained in the range [1.1; 2.1] with only 
4 exception: UGC3851, which has the smallest value 
(as in Paper I); UGC4173, which has a smaller than 
unity value, but still being compatible with the previ- 
ous range at the highest la confidence limit; UGC4325 
and UGC5721 which have larger values. This means that 
the dominant coupling constant for LSBs is the one from 
the gas. 

On the other side, the stellar coupling constant shows 
a two-fold behavior: there galaxies with (3 s tar 3> 1 
(UGC1230, UGC3137, UGC3371, UGC3851, UGC4173, 
UGC4278, UGC8286, UGC9211) and galaxies with 
Pstar < 1 (UGC5721, UGC7524, UGC7603, UGC8837, 
UGC10310). Two systems (UGC1281 and UGC4325) 
have j3 s tar < 1 but their la confidence level are compat- 
ible with values 3> 1. 

Noticeably the systems with the larger j3 sta r are also 
the ones with the smaller stellar mass-to-light ratio. This 
is a warning for a possible degeneracy working among 
these two parameters. Such a degeneracy was already 
discussed in Paper I where we also noticed an anti- 
correlation between the Y* and the only one j3 adopted. 
The fact that we have now broke down this latter in two 
coupling constant ensure us that f3 gas is unaffected by 
any degeneracy, while (3 star is not. This is a major bene- 
fit we have gained by the adoption of the mass component 
decoupling approach; we could add that the main cou- 
pling constant one should rely on is /3 gos . Interestingly 
the gas status is really what it makes the big difference 
between cold and hot dynamical systems which corre- 
spondingly possess cold and hot gas. 

Going more in details of the best fit stellar mass- 
to-light ratios, we have 11 of 15 galaxies compati- 
ble at lcr level with the previously prescribed range; 
one (UGC7524) has an higher value; and three 
(UGC7603,UGC8286, UGC8837) have values smaller 
than 0.5. Among these, the first two were peculiar 
also in Paper I, even if being considered acceptable, 
and UGC7603 in particular still has the smallest value, 
F* = 0.052. Only in two cases (UGC4325 and UGC5721) 
we obtain more compatible values mainly due to their 
large errors. Anyway, in general, in Paper I we obtained 
a better match of F* with the prescribed range in 10 of 
15 cases. 

Finally, the interaction length shows a wider range 
than the clusters galaxies case not being very well con- 
strained as in that case. One particular case being 
UGC3137: even in Paper I it was one of the galaxies 
with the largest interaction length, but now it has a too 
large value, L w 500 kpc, which is difficult to under- 
stand. Anyway, by taking a look to its rotational profile, 
we can see how it satisfies two main requirements: 1. 
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it is the second largest object in the sample (or, bet- 
ter, the second object for which observations of the gas 
component where done up to a larger distance from the 
center); 2. it is the most clear case in which it seems to 
have reached the expected plateau in the rotation curve. 
Lacking more detailed, extended and less disturbed data 
for the galaxies in our sample we cannot conclude any- 
thing; we can only argue that in spiral galaxies, more 
than clusters and the elliptical galaxies, the degeneracy 
among parameters due to limited extension of the data 
could be more important. 

As we did in Paper I, we search for possible scaling 
relations, also for comparing our altern ative approac h 
with more popular theories, like MOND (|Swatersll2010D . 

A correlation among the gas coupling constant (prob- 
ably more related to the global gravitational structure 
properties than the star related parameter) and the ex- 
trapolated central disk surface brightness fJ-o,R (reported 
in Table 3 of Paper I) is possible, see top panel of Fig. [TUJ 

log/3 sas = (-0.32±1.29) + (0.43±Q.96)-log/xo,ii ; (40) 

as long as between the same quantity and the maximum 
rotational velocity, see top panel of Fig. [TT] 

log/3 aas = (-0.30±0.28) + (0.30±0.14) dog V max . (41) 

Excluding only peculiar cases with respect of /3 gas 
(UGC3851, UGC4173, UGC4325, UGC5721) we can see 
that the dispersion around such relations is quite low. 
We also show in the top panel of Fig. [12] the relation 
between (3 gas and the total gas mass of each galaxy: 

log/3 gQS = (0-20 ± 0.06) + (0.04± 0.04)- log M gas ; (42) 

clearly the parameter is quite independent of such quan- 
tity, giving a further confirmation of our preliminary hy- 
pothesis d/3/dr sa 0. In the top panel of Fig.[l3l we show 
the relation with the total baryonic (stars and gas) mass: 

log (3 gas = (0.22 ± 0.06) + (0.03 ± 0.04) • log M bar ; (43) 

We remind that in this case, in the estimation of the stel- 
lar counterpart, we have to take into account the stellar 
mass-to- light ration, Y*, which is actually one of the fit 
parameters. No sensible changes are detected when mov- 
ing from the gas to the total baryonic. 

This does not happen when considering the coupling 
constant with stars; as it is possible that a degeneracy be- 
tween the coupling constant and the stellar mass-to-light 
ratio is working, we consider more properly the quantity 
Pstar ' ^* which appears in the stellar rotational contri- 
bution, finally having (bottom panel of Fig. [TO]) : 

\ogf3 2 star ■ Y, = (-11.69 ± 1.71) + (9.01 ± 1.28) • log/^ij , 

(44) 

and (bottom panel of Fig. [TT]) 

log(3 2 tar -Y, = (1.45± 2.44) + (-0.58 ±1.27)- log V max , 

which holds only for the galaxies with /3 s t ar 3> 1, while 
the others show a more scattered distribution. 

The relation between f3 2 tar ■ Y» and M gas or M bar are 
less evident than the previous case; taking a look to the 
bottom panels of Figs. [T^]- [13] we can only detect (more 
clearly in the M bar case) three sub-groups: two galax- 
ies in the bottom-left corner (/7GG7603 and C/GG8837), 
corresponding to very low values for both j3 s t a r and 



Y*; four galaxies in the middle (C/GG1281, C/GG4325, 
C/GG5721 and J7GG10310), corresponding to low values 
of (3 star] and an almost constant (3 2 tar ■ Y* group made of 
the remaining objects. 

One the other side, the length parameter, L, does not 
show a clear ordered pattern when compared with the 
same quantities as before. 

An even more interesting relation to be tested here 
is the well-know n Baryonic Tully-Fisher (BTF) relation 
(]McGaughll2012[ ). which relates the total baryonic mass 
of spiral galaxies, M bar , with the maximum observed ve- 
locity, Vmax, that in the outer regions becomes approxi- 
mately fiat, Vf and is actually quantified to be: 

M bar =A-Vf, (46) 

with A = 47 ± 6 M© km~ 4 s 4 . In order to make the 
comparison with our results more clear, we can re- write 
the previous in the equivalent expression: 

\ogVj = (1.582 ±0.003) + 0.25 log M bar . (47) 

As it is possible to verify from previously plotted rota- 
tional curves, not all the galaxies in our sample seem to 
have reached the flatness regime in the rotational curve; 
so that the Vf appearing in the previous formulas is more 
exactly V max , the maximum velocity evaluated at the 
maximum distance from the center available from the 
data. 

In our alternative scenario, the total velocity can be 
written as the sum of different terms: 

Vf = V^ star + V Ntgas + V C star + V Cgas , (48) 

where the suffix N and G, ax explained in previous sec- 
tions, are related to the (extended) Newtonian and cor- 
rective terms in our gravitational potential model. 

If we first consider the quantity Vfjheo made up with 
V^^gas entirely derived from direct velocity observations; 
VN,star mixing direct velocity observations and a fit 
parameter, Y* (which also enters in the calculation of 
Mbar); and V^ star and V$ gas completely depending on 
our theoretical model, we obtain: 

logV) = (1.558+ 0.030) + 0.25 log M bar , (49) 

which agree very well with the expected BTF relation. 
We can also check that the BTF relation is verified for 
any of the components of the total velocity (of course, 
we will consider the maximum velocity for any of these). 
For the pseudo-newtonian-term-derived stellar velocity 
we can identify two groups: galaxies with Y* > 2 give 

log V NMar = (0.547 ± 0.067) + 0.25 log M bar , (50) 

and galaxies with Y* < 2 give 

log V NtStar = (1-232 ± 0.043) + 0.25 log M bar . (51) 

The pseudo-newtonian-term-derived gas velocity (ex- 
cluding C/GG4325) gives: 

logV N , gas = (1.132 ± 0.026) + 0.25 log M bar . (52) 

The corrective-term-derived stellar velocity (excluding 
UGC5721, UGC7603, C/GG8837 and C/GG10310) gives: 

\ogV c , gas = (1-092 ±0.064) + 0.25 log M bar . (53) 
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Correlation among the scalar field coupling constants and the central surface brightness of spiral galaxies. 
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The corrective-term-derived gas velocity (excluding 
J7GC3851) gives: 



fogVcstar = (1-448 ±0.051) + 0.25 log M h , 



(54) 



All the previous relations are in Figs. [TH - [15] It is in- 
teresting to note that, for these last two quantities, if 
we relax the condition for which the coefficient of Mf, ar 
has to be equal to 0.25, and we leave it free, we obtain, 
respectively, the values 0.037 and 0.033, which are quite 



consistent with a scenario where the contribution to the 
rotation curve coming from the effective mass produced 
by the modified gravity of the chosen scalar field is prac- 
tically independent of the baryon mass (at least, for the 
family of spiral galaxies considered), so arguing in favor 
of a more universal task for our mechanism. 



5.3. Unified picture 
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Fig. 14. — Baryonic Tully-Fishcr relation for our galaxy sample: contributions from different velocity components. 



In this section we will qualitatively draw some conclu- 
sions from the results that we have derived from the dif- 
ferent gravitational structures discussed above. In par- 
ticular we refer to the new evidence that a scalar field 
can mimic dark matter at different scales by changing its 
properties depending either on the local physical condi- 
tions, namely the matter density, or the matter status. 

As done in Paper I, we consider the relation among 
the scalar field parameters and the gas density/mass of 
any structure in order to check the presence of a global 
correlation among parameters at all scales, from galaxies 
to galaxy clusters. 



For galaxy clusters this is natural to check because the 
gas is the main contribution to the mass (under the as- 
sumption that there is no dark matter, as in the current 
work). The same argument it is not intuitive for galax- 
ies, as gas is not the main component at galactic scale: 
however as we have found that the coupling constant of 
the gas with the scalar field can be larger than the one of 
stars, this component turns out to have a significant im- 
pact on the observational quantities too, sometimes even 
larger than the stellar one. 

Let us start by verifying how the scalar field parame- 
ters are correlated with each other. We remind that one 
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Fig. 15. — Baryonic Tully-Fisher relation for our galaxy sample: total rotational velocity. 



of our main hypothesis is that the coupling constants 
/? are independent of the scale or, more precisely, that 
dj3 1 dr ~ 0. We can verify a posteriori that this hypothe- 
sis is quite well satisfied by taking a look to the top pane 
of Fig. 1161 Here a linear regression of the spiral galaxy 
parameters gives: 

log/3 9as = (0.19 ± 0.05) + (0.04 ± 0.02) • logL ; (55) 
while for clusters of galaxies we obtain: 

log/3 9as = (0.91 ± 0.18) - (0.17 ± 0.06) • logL ; (56) 
and using all data together we finally get: 

log f3 gas = (0.20 ± 0.05) + (0.06 ± 0.02) • log L . (57) 

We can verify how coherently each gravitational struc- 
ture family locate in the parameter space: the clusters of 
galaxies lay in a very well constrained region on the right 
side of top panel in Fig. (THl spiral galaxies are on the op- 
posite side with a larger dispersion mainly due to the 
uncertainties on the interaction length. The only ellipti- 
cal galaxy that we have analyzed is interestingly located 
right in between the spirals and the galaxy cluster regions 
very close to the area where the smallest clusters and/or 
group of galaxies are located. This evidence seems to 
suggest that the properties of the scalar field are related 
to the dynamical properties of the gravitational systems, 
as cold dynamical systems like spirals are clearly sepa- 
rated by hot dynamical systems (elliptical galaxies and 
galaxy clusters). This result is interesting as the con- 
nection between dynamics (especially anisotropy of the 
orbits) and scalar field par ameters has been foun d also 
in other /(r) formulations ( Napolitano et aLll2012[ ). 

For the coupling constant related to stars (in galaxies) 
and galaxies (in clusters), in the bottom panel of Fig. [TBI 
we can see a larger dispersion in the parameter space, 
but also a clear segregation: the left half of the diagram 
is dominated by spiral galaxies, the right half by clusters 
of galaxies, and the central region populated by ellipti- 
cal galaxies (although a larger sample is required) and 
smaller clusters and/or group of galaxies. 

We can eventually reach the same conclusion looking 
at the correlation between scalar field parameters and 
the gas density as in Fig. [T7J when considering the top 
and middle panels of this figure we can see that the pa- 
rameters do not distribute randomly, but seem to follow 



a trend. For the interaction length, the top left corner 
is populated by clusters of galaxies, the bottom right by 
spiral galaxies, and the center possibly by elliptical galax- 
ies and group of galaxies, with a clear linear trend and 
a dispersion which is larger for spiral galaxies than for 
clusters. A linear regression in the log-log space where 
the interaction length L is the dependent variable give, 
for spiral galaxies: 

logL = (10.29 ±4.68) - (1.87 ± 0.96) • log p gas ; (58) 

for clusters of galaxies: 

logL = (5.94 ± 1.37) - (0.78 ±0.35) • logp gas ; (59) 

and for all systems together: 

logL = (9.67 ± 1.04) - (1.74 ± 0.23) • logp 9as 



(60) 

For the gas coupling constant, on the other side, we have 
a linear trend almost compatible with a constant /3 gas , 
thus confirming our previous hypothesis (d/3/dr « 0); in 
particular for spiral galaxies we have: 

log/3 9QS = (-0.31 ±0.57) + (0.11 ± 0.12) -log^ ; (61) 
for clusters of galaxies we have: 

logf3 gas = (-0.04 ±0.41)- (0.11 ± 0.10) -log/Vs ; ( 62 ) 
end for the total sample: 

lo g gas = (1.01 ± 0.18) - (0.16 ± 0.04) • \ogp gas . (63) 

When considering the combination of /3 sta r and Y* or 
f3 g ai, the large spread does not allow to define a clear 
trend and no data regression is produced. 

6. A CONSISTENCY CHECK: STRONG LENSING 
CONVERGENCE MAP 

In previous sections we have extensively tested the dy- 
namical effects of the scalar field theories against the 
stellar dynamics of elliptical galaxies, the gas kinematic 
of spiral systems and the hydrostatical equilibrium of 
galaxy clusters. We have shown that the modified po- 
tential as in Eq. |4] has successfully passed all these 
astrophysical tests by carefully taking care of the cou- 
pling with the different baryonic mass components of 
the gravitating systems (i.e. stars and gas). Most 
of these tests rely, though, on classical mechanics and 
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Fig. 16. — Correlation among the scalar field parameters: spirals arc LSB spiral galaxies; filled squares are clusters of galaxies; filled circle 
is the elliptical galaxy NGC4374. Dashed line is the best fit for clusters only; dotted line is the best fit for spiral galaxies only; dot-dashed 
line is the best fit for the total sample. Objects in brackets are the peculiar cases described in the text and not considered in the fits. 
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Fig. 17. — Correlation among the scalar field parameters and the gas density: spirals are LSB spiral galaxies; filled squares are clusters 
of galaxies; filled circle is the elliptical galaxy NGC4374. Dashed line is the best fit for clusters only; dotted line is the best fit for spiral 
galaxies only; dot-dashed line is the best fit for the total sample. Objects in brackets are the peculiar cases described in the text and not 
considered in the fits. 



hydrostatic equilibrium (i.e. they imply all the same 
kind of physics). One can ask whether the gravita- 
tional effects produced by the scalar field potential can 
be equally effective in reproducing lensing observations, 
which have to be interpreted within general relativity. 
This is a natural pas sage for any alternative or extended 
theories of gravity dBruneton fc Esposito-Faresd 12007 
Angus et all l2006t ~ 



. _ _ - .. .Lubini et al.l 
20111 iBernal et al Jl2011b iBrownstein fc Moffatl l2007f) as 



a failure on accounting for lensing measurements is a too 
strong flaw for a theory to stay stand. 

A complete and exhaustive treatment of gravitational 
lensing predictions and modeling will be presented in a 
forthcoming paper. Here we want to present a simple test 
to visualize whether the proposed approach is suitable to 
extend the scalar field gravity to lensing analyses. 

Due to its intrinsic properties (i.e. the existence of dif- 
ferent coupling constants with different kind of matter) 
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we anticipate that the lensing test might present some 
difficulties from the observational point of view. We re- 
mind that the source of the total potential is given by the 
surface density distribution of baryons (galaxies and gas 
in this case) which in many cluster can be rather noisy 
(e.g. too few cluster members and/or low X-ray emission 
signal-to- noise ratio). As seen for elliptical galaxie, also 
for cluster there are currently not so many systems on 
the market for which it is possible to have high quality 
data from both spectroscopy (i.e., redshift) and X-ray gas 
emissions along with reliable strong and / or weak lensing 
features at the same time. 

As a preliminary test we have performed the strong 
lensing analysis of A383, one of the clusters involved in 
our X-rays analysis as in § 15.11 This cluster is particular 
suitable for our purpose as there are literature measure- 
ments for high resolution data for X-ray gas (as explained 
in the related section here and in Paper I) and recent ob- 
servations from the Cluster Lensing And Supernovae sur- 
vey with Hubble- CLASH project ([Postman et al.ll2011 ) 
which have produced accurate strong dZitrin et al.ll2011l 
Z+ll hereafter) and weak lensing maps (jZitrin et all 
2012). In particular, we have been provided of strong 
lensing convergence (k) map from Z+ll (see their Fig. 4), 
obtained under the standard dark matter paradigm, 
which we have compared against the prediction of the 
scalar field theory. The Z+ll k— map is obtained from 
the multiple-image systems attributed to a galaxy candi- 
date locate at redshift z s — 2.55 (the redshift of A383 is 
zi = 0.189). According to the the standard strong lensing 
equations in Newtonian gravity, the two dimensional sur- 
face mass distribution is rela ted to the lensing-obse rvable 
quantity, the co nvergence k (Schneid er et al.lll9 92). This 
relation writes (|Peacoc kl l2003h : 



n(x,y) 



AkG n DiDu 
c 2 D s 



p(x,y, z)dz 



E 



where 



V(x,y) = / p(x,y,z)dz 



is the two dimensional surface mass density and 

c 2 D„ 



4ttG n DiD 



(64) 



(65) 



(66) 



is the Newtonian critical surface mass density, D s is 
the angular diameter distance to the background source 
galaxy, Di is the angular diameter distance to the lens 
and Dig is the angular diameter distance between the 
lens and the source. All distances are calculated given 
the source and lens redshift, and the fiducial cosmological 
model with Cl m = 0.3, A = 0.7 and h = 0.7. 

In an alternative scenario, we can modify the previous 
formula (jBrownstein fc Moffatil2007tl in: 



K(x,y)= I 4nG } r ^ Dl ^ ls p(x, y, z)dz 



where 



^(x,y)= / G(r)p(x,y,z)dz , 



(67) 
(68) 



is the weighted surface mass density, now not only re- 
lated to the mass distribution, but also to the (alterna- 



tive) gravity theory through the modified gravitational 
constant given in our case by: 



G(r) 
G 



N 



= 1 + 2fi 2 exp(-r/L) 1 + - 



(69) 



Following the same approach of previous sections, we 
have now that the mass density p(x,y, z) is made only of 
the baryonic components (gas and galaxies) , and that the 
scalar field can couple with each of them with different 
strength. 

The K-map from Z+ll is reproduced as light contours 
in our Fig. [TS] where we also overplot the derivation of 
the convergence map from our theoretical model using 
Eqs. [67]- [6l] as black contours. This model, which best 
approximates the original map from Z+ll, corresponds 
to P gas ~ 2.98, P ga i ~ 0.8 and L ~ 590 kpc in Eq. [Ml 

The visual inspection of the map derived in our scalar 
field gravity clearly shows the over-simplification of the 
approach, as the derived k— map is spherically symmet- 
ric while in the Z+ll map the map takes into account 
the local contribution of the massive galaxies in the field 
which perturb the overall convergence of the cluster po- 
tential. This causes a disagreement of the iso— k contours 
in many points, e.g. along the South- West direction , 
where the fairly round behavior of observationally recon- 
structed mass model is partially broken. 

Overall, despite the simple analysis, our model seems 
to recover the main features of the Z+ll map with the 
all the iso— k contours overlapping each other for large 
portions, thus showing a fair correspondence between all 
the plotted mass density level. This is better quantified 
by the 1— dimension k— map in Fig. [18] (right) obtained 
by dividing the surface mass estimations both from data 
and our model in radial bins of 10 kpc. As it can be 
easily seen, the agreement is very good which shows that 
the scalar field potential is compatible also with lensing 
data for A383. 

Importantly, the parameter values which allowed the 
best match to data turned out to be overall compatible 
with the ones derived from the X-ray analysis as in Ta- 
ble[H although there are some slight differences. E.g., the 
coupling with galaxies, f3 ga i, is slightly higher than the 
value derived from the X-ray reconstructed mass anal- 
ysis, but this is not unexpected, as X-ray data are less 
trustworthy in the very inner region of clusters, i.e., for 
R < 70 kpc, where out-of-virial-equilibrium phenomena 
(cooling flow) can alter the results. 

It is also possible that this lensing vs. X-rays discrep- 
ancies can be ascribed to the typical discrepancy found 
also under the dark matter paradigm (see e.g. Comer- 
ford & Natarajan 2007) which might be possibly related 
to intrinsic inconsistency of the technique assumptions 
rather than being physical (e.g. Sereno 2007). 

In summary, we think that this preliminary test reason- 
ably shows that our alternative model can be compatible 
with gravitational lensing. Future steps to improve the 
lensing analysis will include: i) a full theoretical analy- 
sis of gravitational lensing in our alternative scenario; ii) 
the use of a larger sample of galaxy clusters by combin- 
ing X-ray observation and gravitational lensing (where 
possible) in a more self-consistent way; Hi) the inclusion 
of the single galaxies contribution in the same fashion 
of lensing reconstruction models, in order to break the 
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spherical symmetry and to obtain an even better agree- 
ment with data; iv) the extension of the analysis also to 
the weak lensing. 

7. CONCLUSIONS AND DISCUSSION 

In this work we have investigated the dynamical prop- 
erties of several astrophysical systems, from galaxies to 
cluster of galaxies, within the theoretical framework of 
scalar-tensor theories as alternative to the dark mat- 
ter paradigm, and looked for observational signatures in 
galaxy kinematics, X-ray equilibrium and gravitational 
lensing. 

The crucial ingredients of our work are some novel 
properties of the scalar field with respect to previous 
analyses (Paper I) : (i) a new parametrization of the the- 
ory where we have introduced the coupling of the scalar 
field with the different baryon component of the grav- 
itational systems (stars and gas) and (ii) its mass (or 
interaction length) can vary with scale. 

This new parametrization with basically two coupling 
constants, provided allowed the scalar field to mimic dark 
matter, i.e. by reproducing remarkably well observations 
without introducing a new kind of matter and by a mod- 
ification of gravity. 

Going in more details of the results in the different 
astrophysical test we have performed here, results have 
shown that (i) the velocity dispersion of elliptical galax- 
ies can be fit very well by a scalar field, even better than 
assuming a classical Navarro-Frenk- White profile for the 
dark matter component; (ii) a scalar field can fairly well 
reproduce the matter profile in clusters of galaxies, esti- 
mated by X-ray observations and without the need for 
any additional dark matter; (iii) good fits to the rotation 
curves of low surface brightness galaxies are obtained. 
We have finally performed a preliminary test of the mod- 
ified potential derived in the scalar field gravity to repro- 
duce strong lensing convergence map of one of the X-ray 
clusters in our sample. Despite the simplified approach 
adopted, we have obtained a nice match of the k— map 
with the one obtained from a detailed standard DM anal- 
ysis from Z+ll. 

All these results show that scalar field gravity theories 
can be compatible with a wide range of astrophysical 
testa and can be considered as a viable alternative to the 
dark matter. In particular, our results show the possibil- 
ity that the scalar field can couple with ordinary matter 
with different strength (different coupling constants) de- 
pending on the clustering state of matter components in 
the considered gravitational structures and that a pos- 
sible correlation with the evolutionary state of gravita- 
tional systems is possible. 

We point out that this is only the first of a longer se- 
ries of works where we are intended to test our model 



and parametrization by using other different cosmolog- 
ical and astrophysical probes. In future works, on the 
cosmological side, we will study the feasibility of this 
model with the cosmic microwave background radiation 
data and matter power spectrum features; while, on the 
astrophysical side, we will explore the possible corre- 
lations between classical models (NFW or generalized 
NFW profiles) for dark matter and the effective contri- 
bution from our scalar field; for what it concerns spi- 
ral galaxies, it will be interesting to verify if and how 
our model is ab le to reproduce the Universal R otation 
Curve scenario (IPersic et alJH997t ISaluccill20Tl . Even 
more decisive will be to test the model with the gravita- 
tional lensing phenomenon. This is a major probe to rule 
out modified gravity alter natives to dark matter, such 
as no n-relativistic MOND ([Bruneton fc Esposito-Faresd 
120071 ). or to put constraints on other models such as 
TeVeS by requiring the existence of mass ive neutrinos 
(| Angus et all 1200a iFeix et al.1 120081 f2010h . Since our 
model is not a subclass of MOND/TeVeS or the stan- 
dard scalar-tensor theories, it is not ruled out by lens- 
ing observations. In fact, the model studied here, is 
more similar fro m the theoretical point of view, to the 
ones studied in dLubini et all [20Tlt IBernal et at] 120111: 
iBrownstein fc Moffat 2007f). where the authors show that 
all the gravitational predictions obtained in general rela- 
tivity are still valid in alternative theories like chameleon- 
/(i?)-gravities. Due to these reasons we have already 
planned a detailed analysis of strong and weak lensing 
applied to clusters of galaxies, and we believe that the 
possibility of decoupling the different baryonic mass com- 
ponent in a gravitational structure, with different cou- 
pling constant with the chosen scalar field, can make us 
to recover gravitational lensing observations and without 
the need of any massive neutrinos component. 
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TABLE 4 

Clusters of galaxies: Scalar field. Column 1: Name of the cluster. Column 2: Coupling constant of scalar field and galaxy 
component. Column 3: Coupling constant of scalar field and gas component. Columns 4: Scalar field interaction length. 
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TABLE 5 

Spiral galaxies: Scalar field. Column 1: Name of the LSB galaxy. Column 2: Coupling constant of scalar field and star 
component. Column 3: Coupling constant of scalar field and gas component. Columns 4: Scalar field interaction length. 

Column 5: Stellar mass-to-light ratio in the V-band. 
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